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婚 飞 行为 影响 中 华 蜜蜂 性 成 熟 处 女 蜂王 的 基因 表达 


吴 小 波 ， 王 子 龙 , 张 ” 飞 , 石 元 元 ， 曾 志 将 


(江西 农业 大 学 蜜蜂 研究 所 , 南昌 330045 ) 


摘要 : 婚 飞 是 性 成 熟 处 女 蜂 王 与 雄 蜂 交配 过 程 中 的 一 个 重要 前 奏 , 在 该 过 程 中 蜂王 体内 伴随 着 一 系列 重要 的 生理 
变化 。 为 了 探究 中 华 蜜蜂 hpis cerana cerana 处 女 蜂 王 婚 飞 过 程 中 基因 表达 变化 ,本 研究 利用 数字 基因 表达 谱 
( digital gene expression，DGE) 技术 分 析 了 中 华 蜜蜂 性 成 熟 处 女 蜂 王 飞行 与 未 飞行 之 间 的 基因 表达 差异 。 经 DGE 
测序 , 分 别 从 两 个 样品 中 获得 5.98 和 6.01 百 万 条 Clean 标签 。 通 过 分 析 检 测 到 250 个 基因 有 差异 表达 , 其 中 133 
个 基因 在 飞行 蜂王 中 上 调 表达 , 117 个 基因 在 飞行 蜂王 中 下 调 表达 。 这 些 差异 基因 可 以 归 类 到 348 个 功能 性 类 别 
和 142 个 生化 途径 。 结 果 表 明 中 华 蜜蜂 性 成 熟 处 女 蜂 王 在 婚 飞 过 程 中 大 量 基 因 的 表达 发 生 了 变化 。 这 些 结果 为 
进一步 研究 中 华 蜜蜂 蜂王 婚 飞 过 程 中 生理 变化 的 分 子 机 制 提 供 了 重要 的 基因 表达 信息 。 
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Mating flight behaviour affects gene expression in matured virgin queens 


of Apis cerana cerana (Hymenoptera.: Apidae) 

WU Xiao-Bo, WANG Zi-Long, ZHANG Fei, SHI Yuan-Yuan, ZENG Zhi-Jiang”( Honeybee Research 
Institute, Jiangxi Agricultural University, Nanchang 330045, China) 

Abstract: Queen mating flight is an important prelude of sexual matured virgin queens mating with 
drones, which is usually accompanied by a series of important physiological changes in queens ”bodies. 
In order to explore the differences of gene expression following queen mating flight behaviour, we 
analyzed the gene expression differences between the flying matured virgin queens and non-flying matured 
virgin queens of Apis cerana cerana using a high-throughput sequencing method. Through digital gene 
expression (DGE) sequencing, we obtained 5.98 and 6.01 million clean tags from the flying matured 
viregin queen and non-fying matured virgin queen samples, respectively. A total of 230 genes were 
differentially expressed between both, with 133 up-regulated and 117 down-regulated in the flying 
matured virgin queens. These differentially expressed genes can be classified into 348 functional 
categories and involved in 142 biochemical pathways, indicating that there are a large number of genes 
whose expression levels change during mating flight process of matured virgin queens. The results provide 
important gene expression information for further researching the molecular mechanisms of the 
physiological changes during queen mating flight in A. cerana cerana. 

Key words: Apis cerana cerana; queen; flight behaviour; digital gene expression ( DGE ); gene 


expression; mating flight 


蜜蜂 是 与 人 类 以 及 大 自然 关系 密切 的 社会 性 昆 着 分 子 生物 学 技术 的 发 展 , 蜜蜂 分 子 生物 学 研究 也 
忠 , 它 不 仅 可 以 为 人 类 提供 营养 丰 宣 的 蜂 产 品 , 而 取得 了 一 系列 成 果 。 尤 其 是 2006 年 10 月 绘制 完成 
日 通 过 授粉 模式 维持 生态 平衡 ,对 保护 和 维持 生物 西方 蜜蜂 hpis mellifera 基因 组 全 序列 图 谱 后 , 科研 
多 样 性 起 着 重要 作用 。 不 仅 如 此 , 蜜蜂 是 继 果 蝇 后 人员 开 展 了 大 量 的 后 基因 组 学 的 研究 , 对 中 华 蜜 蜂 
又 一 个 重要 的 模式 昆虫 , 被 广泛 应 用 于 学 习 记 忆 、 4. cerana cerana 和 意大利 蜜蜂 4A. mellifera ligustica 
毒 理 安全 以 及 社会 行为 学 研究 ( 兽 志 将 , 2009)。 随 ” 分 子 生物 学 作 了 较为 深入 的 比较 人 研究。 高 通 量 测序 
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技术 (high-throughput sequencing ) 叉 称 “下 一 代 ” 测 
序 技术 , 以 能 一 次 并 行 对 几 十 万 到 几 百 万 条 DNA 
分 子 进行 序列 测定 为 标志 , 可 以 同时 完成 传统 基因 
组 学 研究 (测序 和 注释 ) 以 及 功能 基因 组 学 (基因 表 
达 及 调控 , 基因 功能 , 重 白 /核酸 相互 作用 ) 研 究 
( 刘 掌 襄 , 2012; 视 华 良 等 , 2012)。 数 字 基 因 表 达 
谱 ( digital gene expression，DGE ) 就 是 利用 新 一 代 高 
通 量 测序 技术 和 高 性 能 计算 分 析 技 术 来 全 面 、 快速 
地 检测 某 一 物种 特定 组 织 特定 条 件 下 的 基因 表达 情 
况 , 该 技术 目前 已 成 功 应 用 于 蜜蜂 分 子 生物 学 研究 
当中 , 并 发 现 不 同 蜂 种 同一 生理 时 期 以 及 同一 蜂 种 
不 同 生 理 时 期 蜜蜂 的 基因 表达 存在 差异 ( Kocher et 
al., 2008, 2010; Lu ef al., 2011; Wang et al., 
2012 ) 。 

蜂王 是 蜜蜂 蜂 群 当中 的 主要 繁殖 个 体 , 其 生理 
情况 随 日 龄 和 交配 行为 而 改变 。Kocher 等 (2008， 
2010 ) 比较 分 析 了 处 女 蜂王 、 刚 授精 蜂王 以 及 产 卵 
蜂王 的 转录 组 , 分 析 发 现 它们 之 间 存 在 大 量 差 异 ， 
而 且 处 女 蜂 王 授 精 方式 对 其 基因 表达 也 有 一 定 的 影 
啊 。 然 而 , 目前 对 于 中 华 蜜蜂 处 女 蜂 王 的 分 子 生物 
学 比较 研究 报道 较 少 , 引 人 关 注 的 是 , 性 成 熟 处 女 
蜂王 在 蜂 群 内 不 会 吸引 雄 蜂 进行 交配 , 但 在 空中 进 
行 婚 飞 时 却 能 吸引 雄 蜂 与 之 交配 , 表明 处 女 蜂王 在 
婚 发 时 可 能 产生 了 某 些 内 在 生理 变化 从 而 吸引 了 雄 
蜂 的 进行 交配 。 本 人 研究 利用 数字 基因 表达 谱 技 术 比 
较 分 析 婚 飞 对 中 华 蜜蜂 性 成 谎 期 处 女 蜂 王 mRNA 
表达 的 影响 , 结果 如 下 。 


1 材料 与 万 法 


1.1 供 试 昆虫 

试验 蜂 群 为 江西 农业 大 学 蜜蜂 研究 所 饲养 的 中 华 
蜜蜂 4，cerana cerana。 群 势 均 为 5 脾 蜂 , 蜂 脾 相称 。 
1.2 样品 的 采集 

按照 标准 的 人 工 育 王 方法 ( 曾 志 将 , 2009 ) 培育 
中 华 蜜蜂 蜂王 , 待 蜂王 性 成 熟 后 进行 单 雄 人 工 授 
精 。 并 以 此 蜂 群 为 母 本 和 群 提 供 小 幼虫 , 进一步 进行 
人 工 育 王 , 将 新 出 房 的 处 女王 内 入 王 笼 并 放 人 无 王 
区 。 处 女 蜂 王 出 房 12 d 后 , 将 其 中 6 头 处 女 蜂王 分 
别 放 入 含有 几 百 涉 雄 蜂 和 工蜂 的 网 单 ( 长 x 宽 x 融 
=3mx4 mx5 m) 中 ,飞行 20 min 后 采集 飞行 处 
女 蜂王 , 立即 与 为 外 6 头 同 批 次 未 飞行 的 处 女 蜂 王 
一 起 放 入 -80Y 冷冻 处 理 5 min,， 用 解剖 剪 从 蜂王 
侧面 解 齐 并 剔除 蜂王 体内 的 消化 道 , 每 2 头 蜂 王 为 


一 个 样品 , 放 入 - 80 冰箱 保存 , 待 提取 RNA。 
1.3 创建 DGE 文库 并 测序 

用 SV 总 RNA 解 离 系统 (Promega 公司 , 美国 ) 
分 别提 取 中 华 蜜蜂 性 成 熟 处 女 蜂 王 每 个 样品 的 总 
RNA, 并 进行 RNA 质量 检测 。 竺 样品 检验 合格 后 ， 
分 别 将 中 华 蜜蜂 性 成 熟 处 女 蜂 王 的 总 RNA 平行 愉 
品 混合 , 制备 飞行 蜂王 和 未 飞行 蜂王 试验 样品 。 分 
别 取 6 pg 总 RNA, 利用 Oligo( dT) 磁 珠 定 集 总 RNA 
中 的 mRNA, 并 以 Oligo(dT) 为 引导 逆转 录 为 双 链 
cDNA , 采用 4 碱 基 识别 酶 We 亚 酶 切 双 链 cDNA， 
连接 上 JUlumina adaptor 1， 利 用 MneTI 酶 切 3 端 
CATG 下 游 17 bp 碱 基 ,并 在 3“ 端 连接 上 Humina 
adaptor 2。 再 加 入 Primer GX1 和 Primer GX2 进行 
15 个 循环 的 PCR 扩 增 。 扩 增产 物 通过 6% TBE 
PAGE 胶 回 收 95 碱 基 条 带 ， 纯 化 后 的 族 段 用 
Ilumina HiSeq ”2000 测序 仪 (Jlumina 公司 , 美国 ) 
进行 测序 。 测 序 产生 数 百 万 条 原始 Reads，Reads 
的 测序 读 长 为 49 bp。 
1.4 DGE 标签 与 意大利 蜜蜂 基因 序列 的 比 对 分 析 

测序 得 到 的 原始 图 像 数据 经 读 取 测 序 峰 图 转化 
为 序列 数据 , 称 之 为 原始 数据 或 原始 reads， 以 fastq 
文件 格式 存储 , 里 面 存储 reads 的 序列 以 及 reads 的 
测序 质量 。 为 了 获得 理想 的 序列 数据 ,原始 数据 经 
过 去 除 空 载 reads( 只 含 3' 接 头 而 找 不 到 标签 序列 的 
reads) 、 去 除 3 "接头 序列 (由 于 标签 只 有 21 nt 而 测 
序 read 长 度 为 49 nt, 原始 read 上 市 有 一 段 3' 接 头 
序列 )、 去 除 低 质量 标签 (含有 未 知 碱 基 N 的 标 
签 )、 去 除 长 度 过 小 或 过 大 的 标签 , 取 长 度 为 21 nt 
的 标签 、 去 除 拷贝 数 为 1 的 标签 (可 能 是 测序 错 
误 ) ， 从 而 获得 包含 CATG 的 21 bp 标签 序列 的 
Clean 标签 。 

从 NCBI (ftp://ftp. necbi. nih. gov/genomes/ Apis_ 
mellifera ) 下 载 意 大 利 蜜蜂 所 有 基因 的 mRNA 序列 ， 
利用 软件 检索 mRNA 上 所 有 的 CATG 位 点 ,生成 
CATG +17 nt 碱 基 的 参考 标签 数据 库 。 然 后 将 试验 
获取 的 全 部 Clean Tag 与 参考 标签 数据 库 进行 比 对 ， 
允许 最 多 一 个 碱 基 错 配 , 对 其 中 唯一 比 对 到 一 个 基 
的 标签 (unambiguous tags) 进行 基因 注释 , 统计 每 
个 基因 对 应 的 原始 Clean Tag 数 , 然后 对 原始 Clean 
Tag 数 做 标准 化 (TPM ) 处 理 , 标准 化 方法 为 : 每 个 
基因 包含 的 原始 Clean Tag 数 /该 样本 中 总 Clean 
Tag 数 x1 000 000, 从 而 获得 标准 化 的 基因 表达 量 。 
1.5 ”差异 表达 基因 的 筛选 

根据 Audic 和 Claverie (1997 ) 描述 的 数字 化 基 
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因 表达 谱 差 异 基因 检测 方法 ,筛选 飞行 蜂王 和 未 飞 
行 蜂王 两 样本 间 的 差异 表达 基因 。 在 多 重 假 设 检验 
校正 和 分 析 中 FDR 被 用 来 确定 尸 值 的 域 值 。 我 们 
用 FDR <0.001 且 倍 数 差 异 在 2 倍 及 以 上 作为 临界 
值 去 分 析 差 异 表达 基因 。 在 此 基础 上 进行 Cene 
Ontology 功能 显著 性 富 集 分 析 和 KEGG Pathway 显 
闭 性 富 集 分 析 。 
1.6 数据 统计 与 分 析 

参照 Audic 和 Claverie (1997 ) 描述 的 方法 , 假 
设 观 测 到 茶 基 因 对 应 的 Clean 标签 数 为 x*, 已 知 在 
文库 中 , 每 个 基因 的 表达 量 只 占 所 有 基因 表达 量 的 
一 小 部 分 , 在 此 情况 下 , p(x) 的 分 布 服 从 沾 松 分 布 : 


p(z) = 和。 (》 为 基因 A 的 真实 转录 数 ) 


已 知 , 样本 1 总 Clean 标签 数 为 Nl, 样本 2 总 
Clean 标签 数 为 N2, 某 基 因 A 在 样本 1 中 对 应 的 
Clean 标签 数 为 x, 在 样本 2 中 对 应 的 Clean 标签 数 
为 y, 则 此 基因 A 在 两 样本 中 表达 量 相 等 的 概率 由 
以 下 公式 来 计算 : 


2 pCi % ) 





或 2 x (1 - pi 6) ) (如 果 王 pG x) > 0.5) 


N +y)! 
p(y1 4) = (RF) 
bd bd AN 


假设 挑选 了 M 个 差异 表达 基因 ,其 中 0 个 是 
真正 有 差异 表达 的 基因 , 另外 N 个 是 其 实 没 有 差异 
表达 的 基因 , 为 假 阳 性 。 错 误 比 例 S = N/M 平均 而 
言 不 应 超过 某 个 临界 值 , 在 统计 时 预先 设 定 其 FDR 
值 不 能 超过 1% 。 在 获得 差异 检验 的 FDR 值 时 , 也 
根据 基因 的 表达 量 (RPKM 值 ) 计算 此 基因 在 不 同 
样本 间 的 差异 表达 倍数 。 差 异 倍数 越 大 ，FDR 值 越 
小 ,就 表明 表达 差异 起 显著 。 在 此 分 析 中 ， 差异 表 
达 基 因 被 定义 为 FEDR 三 1%o 上 且 倍 数 差 异 在 2 倍 以 上 
( 合 2 倍 ) 的 基因 。 


2 结果 与 分 析 


2.1 中 华 蜜蜂 飞行 蜂王 与 未 飞行 蜂王 DGE 序列 
分 析 

分 别 构建 中 华 蜜蜂 飞行 蜂王 (flying matured 
virgin queen,，FMVQ、Q ) 和 未 飞行 蜂王 (non-flying 
matured virgin queen,，NFMVQ ) 的 DGE 数据 库 并 测 
序 , 分 别 获 得 6 104 196 和 6 123 579 条 原始 标签 。 


过 滤 排 除 低 质量 的 标签 后 ， 每 个 数据 库 获 得 的 总 
Clean 标签 数 分别 是 5 980 022 和 6 010 094( 表 1)， 
这 些 Clean 标签 分 别 占 总 标签 数 的 97. 97% 和 
98.15% (图 1)。 找 贝 数 大 于 100 的 标签 的 总 数 在 
Clean 标签 之 中 所 占 的 比例 超过 77% , 但 是 它们 的 
种 类 在 Clean 标签 种 类 中 的 分 布 不 超过 6% , 拷贝 
数 在 2 ~5 之 间 的 标签 种 类 在 Clean 标签 种 类 中 的 
分 布 比 较 广 泛 ( 图 2)。 表 达 谱 测序 饱和 度 分 析 用 于 
检查 检测 到 的 基因 数目 是 否 随 着 测序 量 ( 标 签 数 
量 ，total tag number) 的 增加 而 增加 并 且 赵 于 饱和 。 
如 图 3 所 示 ， 当 两 数据 库 的 测序 量 达 到 接近 2 M 
时 , 检测 到 的 基因 数目 趋 于 饱和 。 
2.2 DGE 序列 比 对 到 意大利 蜜蜂 基因 序列 

由 于 目前 GenBank 中 仪 有 意大利 蜜蜂 的 基因 
组 和 基因 序列 , 而 中 华 蜜 蜂 的 基因 组 序列 数据 尚未 
完成 。 而 且 中 华 蜜蜂 和 意大利 蜜蜂 在 蜜蜂 属 中 是 亲 
缘 关 系 较 近 的 两 个 种 , 理论 上 两 者 间 基 因 组 序列 差 
异 应 该 很 小 , 因此 , 在 本 人 研究 中 我 们 以 意大利 蜜蜂 
的 基因 序列 作为 参考 序列 来 分 析 基 因 的 表达 情况 。 
在 Clean 标签 中 , 两 个 样品 中 能 够 比 对 到 参考 意 大 
利 蜜蜂 基因 序列 的 标签 总 数 分 别 为 1 896 575 和 
1 857 672, 分 别 占 31. 72% 和 30.91% 。 能 比 对 到 
参考 基因 的 标签 的 种 类 数 及 百分比 分 别 是 44 318 
( 占 37.18%) 和 45 775( 占 39.37%)。 其 中 ,所 有 
Clean 标签 种 类 的 36. 51% 和 38.66% 能 比 对 到 单 
一 序列 (图 4)。 
2.3 “飞行 蜂王 与 未 飞行 蜂王 间 差 异 表达 基因 

为 了 检测 飞行 蜂王 与 未 飞行 性 成 熟 处 女 蜂王 间 
的 差异 表达 基因 , 我 们 根据 Audic 和 Claverie 
(1997 ) 描述 的 数字 化 基因 表达 谱 差 异 基因 检测 方 
法 ,筛选 两 样本 间 的 差异 表达 基因 。 在 两 组 蜂王 数 
据 库 中 , 总 共 检 测 到 250 个 差异 表达 基因 (图 5)， 
飞行 蜂王 与 不 飞行 蜂王 相 比 有 133 个 上 调 表达 基因 
和 117 个 下 调 表 达 基 因 。 其 中 , 类别 最 多 的 为 核糖 
体 蛋 白 、 线粒体 核糖 体 蛋 日 、 细 胞 色素 以 及 丝氨酸 
和 蛋白酶 。 男 外 还 发 现 : 在 飞行 的 中 华 蜜蜂 处 女 蜂 王 
中 气味 结合 蛋 日 9( 0BP9 ) 基因 的 表达 量 比 未 飞行 
的 处 女 蜂 王 要 高 ( 表 2)。 
2.4 差异 基因 的 GO 功能 分 析 与 Pathway 富 集 
分 析 

用 G0 分 类 分 析 这 些 差异 表达 基因 上 所 属 的 功能 
类 别 。 对 于 250 个 差异 表达 基因 , 仪 有 117 个 基因 
有 GO ID, 并 按照 细胞 成 分 、 分 子 功能 和 生化 过 程 
3 个 大 类 , 共 分 为 348 个 功能 性 类 别 。 当 把 这 些 差 
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异 表 达 基 因 与 CO 数据 库 进 行 比 对 并 与 整个 基因 组 。 显著 宣 集 (P<0.05)。 
进行 比较 时 , 发 现 作用 于 糖 基 键 的 水 解 酶 相关 基因 


表 1 中 华 蜜蜂 飞行 蜂王 (FMVQ) 与 未 飞行 蜂王 (NFMVQ ) 数字 化 基因 表达 谱 序 列 分 析 
Table 1 The DGE sequencing in the flying matured virgin queen (了 了 MVYQ ) and non-flying matured virgin queen 


(NFMVEQ) of Apis cerana cerana 


NFMVQ FMVQ 
类 别 - F rp 
ntl 标签 总 数 标签 种 类 标签 总 数 标签 种 类 
Total tag number Distinct tag number Total tag number Distinct tag number 
原 
其 数据 6 123 579 228 005 6 104 196 241 625 
Raw data 
Cl 未 签 
on 6 010 094 116 275 5 980 022 119 213 
Clean tag 


比 对 到 参考 基因 的 所 有 标签 


All tags mapping to the reference genes 


y = 不 
Wan 1 813 329 (30.17% ) 44 947 (38.66%) 1 850 025 (30.94% ) 43 522 (36.51% ) 


Unambiguous tags mapping to the unique reference gene 


比 对 到 基因 组 
Mapping to genome 
未 比 对 上 的 标签 


Unknown tags 


表 中 括号 内 的 数据 为 其 对 应 的 标签 总 数 或 标签 种 类 占 Clean 标签 的 比例 。The numbers in parentheses indicate the percentage of total tag number or 


distinct tag number in clean tags. 


1 857 672 (30.91% ) 45775 (39.37%) 1 896 575 (31.72% ) 44 318 (37.18% ) 


1 228 394 (20.44%) 32 647 (28.08%) 1056 523 (17.67%) 32 858 (27.56% ) 


2 924 028 (48.65% ) 37 853 (32.55% ) 3 026 924 (50.62%) 42 037 (35.26% ) 


标签 分 类 Tag classification 标签 分 类 Tag classification 







加 含 N 的 标签 数 Tags containing N (2 406, 1.00%) 
国 仅 有 接头 Only adaptors (25, 0.01%) 

天 拷贝 数 <2 Copy number<2 (119 981, 49.66%) 
转 干净 标签 Clean tags (119 213, 49.34%) 







加 含 N 的 标签 数 Tags containing N (4 163, 0.07%) 
国 仅 有 接头 Only adaptors (30, 0.00%) 

着 拷贝 数 <2 Copy number<2 (119 981, 1.97%) 
国 干净 标签 Clean tags (5 980 022, 97.97%) 


标签 分 类 Tag classification 标签 分 类 Tag classification 







图 含 N 的 标签 数 Tags containing N (4 013, 0.07%) 
图 仅 有 接头 Only adaptors (32, 0.00%) 

国 拷贝 数 <2 Copy number<2 (109 440, 1.79%) 

转 干净 标签 Clean tags (6 010 094, 98.15%) 


因 含 N 的 标签 数 Tags containing N (2 264, 0.99%) 
图 仅 有 接头 Only adaptors (26, 0.01%) 

着 拷贝 数 <2 Copy number<2 (109 440, 48.00%) 
转 干净 标签 Clean tags (116 275, 51.00%) 





图 1 中 华 蜜蜂 飞行 蜂王 (FMVQ) 与 未 飞行 蜂王 (NFMVQ) 中 不 同类 别 标签 在 原始 标签 中 的 分 布 
Fig. 1 Distribution of raw tags over different tag abundance categories in the flying matured virgin queen (FMVOQ) 
and non-flying matured virgin queen (NFMVQ) of Apis cerana cerana 
A: 飞行 蜂王 标签 总 数 分 布 Distribution of total tags in FMVQ; B: 飞行 蜂王 标签 种 类 的 分 布 Distribution of distinct tags in FMVQ; C: 未 飞行 蜂王 
标签 总 数 分 布 Distribution of total tags ip NFMVQ; D: 未 飞行 蜂王 标签 种 类 的 分 布 Distribution of distinct tags in NFMVQ. 括号 内 的 数字 表示 每 个 
类 型 的 标签 的 总 个 数 及 所 占 的 百分比 。The numbers in parentheses indicate the total numbers and the percentage of tags of each type among the total 


raw tags. 















标签 拷贝 数 Tag copy number 


[2, 5] (181 092, 3.03%) 

[6, 10] (144 002, 2.41%) 
[11, 20] (203 657, 3.41%) 
[21, 50] (395 590, 6.62%) 
[51, 100] (405 212, 6.78%) 
[>100] (4 650 469, 77.77%) 


标签 拷贝 数 Tag copy number 


国 [2, 5] (174 327, 2.90%) 

国 [6, 10] (142 319, 2.37%) 
国 [11, 20] (205 273, 3.42%) 
国 [21, 50] (399 107, 6.64%) 
国 [51, 100] (407 501, 6.78%) 
国 [>100] (4 681 567, 77.90%) 
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B 


标签 拷贝 数 Tag copy number 







国 [2, 5] (62 125, 52.11%) 
国 [6, 10] (18 921, 15.87%) 
[11, 20] (13 915, 11.67%) 
国 [21, 50] (12 333, 10.35%) 
国 [51, 100] (5 710, 4.79%) 
国 [>100] (6 209, 5.21%) 





这 


签 拷 贝 数 Tag copy number 


[2, 5] (59 425, 51.11%) 
[6, 10] (18 673, 16.06%) 
[11, 20] (13 997, 12.04%) 
[21, 50] (12 445, 10.70%) 
[51, 100] (5 783, 4.97%) 
[>100] (5 952, 5.12%) 





图 2 中华 蜜蜂 飞行 蜂王 和 未 飞行 蜂王 Clean 标签 拷贝 数 分 布 统计 
Fig. 2 ”Distrbution of clean tag copy number in the flying matured virgin queen (FMVQ) and 


non-flying matured virgin queen (NFMVQ) of Apis cerana cerana 

A: 飞行 蜂王 Clean 标签 总 数 分 布 Distribution of total clean tags in FMVQ; B: 飞行 蜂王 Clean 标签 种 类 的 分 布 Distribution of distinct clean tags in 
FMVQ; C: 未 飞行 蜂王 Clean 标签 总 数 分 布 Distribution of total clean tags in NFMVQ; D: 未 飞行 蜂王 Clean 标签 种 类 的 分 布 Distribution of distinet 
clean tags in NFMVQ. 在 方 括号 内 的 数字 表示 某 一 特定 种 类 标签 的 拷贝 数 的 范围 ; 在 括 弧 内 的 数字 表示 在 这 一 类 别 中 标签 的 总 个 数 及 占 的 比 


例 。Numbers in the square brackets indicate the range of copy numbers for a specific category of tags ，while the numbers in the parentheses show the total 


numbers and percentage of tags in that category. 


人 


基因 百分比 
Percentage of genes 


100 


80 


干净 标签 匹配 上 的 基因 
Genes mapped by all clean tags 


明确 干净 标签 匹配 上 的 基因 


Genes mapped by unambiguous clean tags 





10 20 30 40 


总 标签 数 (x100 K) 


Total tag number (x100 K) 


B 
100 一 干净 标签 匹配 上 的 基因 
Genes mapped by all clean tags 
---- 明确 干净 标签 匹配 上 的 基因 
80 Genes mapped by unambiguous clean tags 


CN 
(em 


基因 百分比 
Percentage of genes 

上 

S 


20 





50 60 0 10 20 30 40 50 60 


Total tag number (x100 K) 


图 3 中 华 蜜蜂 飞行 蜂王 ( A) 与 未 飞行 蜂王 (B) Clean 标签 测序 饱和 度 分 析 
Fig. 3 Sequencing saturation analysis of clean tags in the flying matured virgin queen (FMVQ) (A) and 


non-flying matured virgin queen (NFMVQ) (B) of Apis cerana cerana 


A: 飞 行 蜂王 Clean 标签 总 数 比 对 到 意大利 蜜蜂 基因 和 基因 组 序列 的 分 布 图 Distribution of total clean tags in FMVQ mapping to genes and genome of 
Apis mellifera ligustica; B: 飞行 蜂王 Clean 标签 种 类 比 对 到 意大利 蜜蜂 基因 和 基因 组 序列 的 分 布 图 Distribution of distinct clean tags in FMVQ 
mapping to genes and genome of A. mellifera ligustica; C: 未 飞行 蜂王 Clean 标签 总 数 比 对 到 意大利 蜜蜂 基因 和 基因 组 序列 的 分 布 图 Distribution 
of total clean tags in NFMVQ mapping to genes and genome of A. mellifera ligustica; D: 未 飞行 蜂王 Clean 标签 种 类 比 对 到 意大利 蜜蜂 基因 和 基因 组 
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标签 分 类 Tag classification 


(PM:Perfect Match) 
(MM:Miss Match) 


标签 分 类 Tag classification 


(PM:Perfect Match) 
(MM:Miss Match) 


PM(Sense) 1 tag->1 gene (1 186 953, 19.85%) 

加 PM(Sense) 1 tag->n gene (22 236, 0.37%) 

二 1MM(Sense) 1 tag->1 gene (442 463, 7.40%) 

国 1 MM(Sense) 1 tag->n gene (20 809, 0.35%) 

国 PM(AntiSense) 1 tag->1 gene (158 652, 2.65%) 
PM(AntiSense) 1 tag->n gene (1747, 0.03%) 

1 MM(AntiSense) 1 tag->1 gene (61 957, 1.04%) 

1 MM(AntiSense) 1 tag->n gene (1 758, 0.03%) 
PM Genome 1 tag->1 position (357 279, 5.97%) 
PM Genome 1 tag->n position (8 826, 0.15%) 

1 MM Genome 1 tag->1 position (354 517, 5.93%) 
1 MM Genome 1 tag->n position (335 901, 5.62%) 
Unknown Tag (3 026 924, 50.62%) 


PM(Sense) 1 tag->1 gene (1 163 348, 19.36%) 

图 PM(Sense) 1 tag->n gene (23 780, 0.40%) 

二 1MM(Sense) 1 tag->1 gene (436 251, 7.26%) 

国 1 MM(Sense) 1 tag->n gene (16 801, 0.28%) 
PM(AntiSense) 1 tag->1 gene (150 753, 2.51%) 
PM(AntiSense) 1 tag->n gene (2 100, 0.03%) 

1 MM(AntiSense) 1 tag->1 gene (62 977, 1.05%) 
1 MM(AntiSense) 1 tag->n gene (1 662, 0.03%) 
PM Genome 1 tag->1 position (365 159, 6.08%) 
PM Genome 1 tag->n position (8 086, 0.13%) 

1 MM Genome 1 tag->1 position (357 852, 5.95%) 
1 MM Genome 1 tag->n position (497 297, 8.27%) 
Unknown Tag (2 924 028, 48.65%) 








标签 分 类 Tag classification 


PM(Sense) 1 tag->1 gene (20 334, 17.06%) 

团 PM(Sense) 1 tag->n gene (274, 0.23%) 

1 MM(Sense) 1 tag->1 gene (8 770, 7.36%) 

国 1 MM(Sense) 1 tag->n gene (287, 0.24%) 
PM(AntiSense) 1 tag->1 gene (10 484, 8.79%) 
PM(AntiSense) 1 tag->n gene (130, 0.11%) 

1 MM(AntiSense) 1 tag->1 gene (3 934, 3.30%) 
1 MM(AntiSense) 1 tag->n gene (105, 0.09%) 
PM Genome 1 tag->1 position (15 292, 12.83%) 
PM Genome 1 tag->n position (427, 0.36%) 

1 MM Genome 1 tag->1 position (14 857, 12.46%) 
1 MM Genome 1 tag->n position (2 282, 1.91%) 
Unknown Tag (42 037, 35.26%) 

(PM:Perfect Match) 

(MM:Miss Match) 


标签 分 类 Tag classification 


略 PM(Sense) 1 tag->1 gene (21 621, 18.59%) 
国 PM(Sense) 1 tag->n gene (290, 0.25%) 
1 MM(Sense) 1 tag->1 gene (8 955, 7.70%) 
国 1 MM(Sense) 1 tag->n gene (286, 0.25%) 
国 PM(AntiSense) 1 tag->1 gene (10 478, 9.01%) 
PM(AntiSense) 1 tag->n gene (136, 0.12%) 
1 MM(AntiSense) 1 tag->1 gene (3 893, 3.35%) 
1 MM(AntiSense) 1 tag->n gene (116, 0.10%) 
PM Genome 1 tag->1 position (15 455, 13.29%) 
PM Genome 1 tag->n position (435, 0.37%) 
1 MM Genome 1 tag->1 position (14 618, 12.57%) 
1 MM Genome 1 tag->n position (2 139, 1.84%) 
国 Unknown Tag (37 853, 32.55%) 

(PM:Perfect Match) 

(MM:Miss Match) 


[1 | II 


图 4 ” 中华 蜜蜂 飞行 蜂王 与 未 飞行 蜂王 Clean tag 比 对 统计 


Fig. 4 Distribution of stat of clean tag alignment in the flying matured virgin queen (FMVQ) and 


non-flying matured virgin queen (NFMVQ) of Apis cerana cerana 


序列 的 分 布 图 Distribution of distinct clean tags in NFMVQ mapping to genes and genome of A. mellifera ligustica. 


表 2 中 华 蜜蜂 飞行 蜂王 与 未 飞行 蜂王 差异 表达 基因 
Table 2 The differentially expressed genes in the flying matured virgin queen (FMVQ) and 


non-flying matured virgin queen (NEMVQ ) of Apis cerana cerana 


基因 

Gene 
核糖 体 蛋白 S13 基因 RpS13 gene 
核糖 体 蛋白 S26 基因 RpS26 gene 
核糖 体 蛋 白 L19 基因 RpL19 gene 
核糖 体 蛋 白 L23A 基因 RpL23A gene 
线粒体 核糖 体 蛋 白 S17 基因 mRpS17 gene 
线粒体 核糖 体 蛋 白 L16 基因 mRpL16 gene 
细胞 色素 P450 9Q1 基因 CYP9Q1 gene 
细胞 色素 P450 305D1 基因 CYP305D1 gene 
细胞 色素 P450 6AS10 基因 CYP6AS10 gene 
丝氨酸 蛋白 酶 22 基因 SP22 gene 
丝氨酸 和 蛋白酶 35 基因 SP35 gene 
丝氨酸 蛋白 酶 40 基因 SP40 gene 
气味 结合 蛋白 9 基因 OBP9 gene 
丝氨酸 蛋白 酶 23 基因 SP23 gene 


+ : 基因 表达 上 调 Up-regulated; - : 基因 表达 下 调 Down-regulated. 


生成 强度 Raw intensity 调控 模式 
FMVQ NFMVQ Regulation model 
4 S01 2 070 十 
1 $529 755 十 
568 184 
2 391 1 139 十 
92 43 
197 83 十 
260 125 十 
244 118 + 
166 8 十 
2 971 之 71 十 
44 3 十 
96 25 十 
67 11 十 
46 175 = 
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表 3 KEGG Pathway 差异 显著 性 富 集 分 析 结 果 
Table 3 Result of KEGG Pathway enrichment analysis 
注释 到 该 通路 的 差异 注释 到 该 通路 的 
通路 表达 基因 的 数目 所 有 基因 的 数目 P 值 Q@ 值 通路 ID 号 
Pathway Number of DEGs with Number of genes with P-value Q-value Pathway ID 
pathway annotation pathway annotation 
维生素 消化 上 
WU 8 56 4.03E -05 0. 005716021 ko04977 
Vitamin digestion and absorption 
胸腺 分 小 12 148 0.000161257 0.011449233 ko04972 
Pancreatic secretion 
化 上 
朋 肪 消化 吸收 7 57 0.000326153 0.015437928 ko04975 
Fat digestion and absorption 
旧 测 人 滑 ¥ 62 0.000550335 0.019536907 ko00561 
Glycerolipid metabolism 
代谢 通路 41 1100 0. 001110727 0.031544647 ko01100 


Metabolic pathways 


4 国 上 调 基 因 Up-regulated genes 
国 下 调 基 因 Down-regulated genes _ 
国 非 差 异 表达 基因 Non-differentially expressed genes ， 


3 - 


et 
- 


Relative expression level of genes in FMVQ 


飞行 蜂王 基因 相对 表达 水 平 4g FMVQ<TPM>) 
a 


-2 





一 | 0 1 2 3 二 


未 飞行 蜂王 基因 相对 表达 水 平 4g NFMVQ<TPM>) 
Relative expression level of genes In NFMVQ 


图 5 中 华 蜜蜂 飞行 蜂王 与 未 飞行 蜂王 基因 表达 水 平 
Fig. S$ The expression level of genes in the flying matured 


virgin queen (FMVQ) and non-flying matured virgin queen 
(NFMVQ) of Apis cerana cerana 
图 中 TPM (transcripts per million clean tags) 是 一 个 标准 化 的 指标 , 指 
每 100 万 Clean tag 中 包含 该 种 转录 本 的 拷贝 数 。 通 过 立 值 FDR < 
0.001 和 |log, Ratiol > =1 (Ratio: 飞行 蜂王 /未 飞行 蜂王 ) 来 确定 有 
表达 差异 的 基因 ; 图 中 所 示 的 红 点 则 表示 在 飞行 蜂王 中 上 调 的 基 
因 , 绿色 点 则 表示 在 飞行 蜂王 中 下 调 基因 ， 蓝 色 则 表示 在 两 种 蜂王 
中 表达 差异 不 显著 的 基因 。Transcripts per million clean tags (TPM ) 
is a normalized index which denotes the number of a certain tag per 
million clean tags. The false-discovery rate (FDR) of <0.001 and 
|log, Ratio| 三 1 (Ratio: FMVQ/NFMVQ) are used as the criterion to 














identify the differentially expressed genes. The red shows the up-regulated 
genes in FMVQOQ, the green represents the down-regulated genes in 
FMVQ, and the blue shows the non-differentially expressed genes 
between NFMVOQ and FMVQ. 


进一步 研究 这 些 差异 表达 基因 所 涉及 的 生化 调 
节 途 径 , 我 们 把 这 些 差异 表达 基因 与 KEGG 数据 库 
进行 比 对 , 发 现 这 些 基因 能 归 类 到 142 个 生化 途 
径 , 当 将 这 些 差异 基因 与 整个 基因 组 背景 进行 比较 
时 , 发 现 与 维生素 消化 与 吸收 , 胰腺 分 泌 , 脂肪 消 
化 与 吸收 ,甘油 代谢 和 代谢 通路 有 关 的 基因 显著 性 
镶 集 (Q <0.05)。 


3 讨论 


DGE 技术 是 利用 新 一 代 高 通 量 测序 技术 和 高 
性 能 计算 分 析 技 术 , 能 够 全 面 、 经 济 、 快速 地 检测 
某 一 物种 特定 组 织 在 特定 状态 下 相关 基因 表达 情 
况 , 目前 已 广泛 应 用 于 基础 研究 , 该 技术 也 已 经 被 
运用 于 蜜蜂 生物 学 特性 的 分 子 机 制 研 究 。Liu 等 
(2011) 利 用 DGE 分 析 了 意大利 蜜蜂 哺育 蜂 和 采集 
蜂 差 异 表 达 基 因 。Wang 等 (2012 ) 研究 了 中 华 蜜蜂 
蜂王 和 工蜂 之 间 的 差异 表达 基因 。Shi 等 (2012 ) 研 
究 了 营养 杂交 对 基因 表达 的 影响 ,等 等 。 本 研究 利 
用 DGE 分 析 了 中 华 蜜蜂 飞行 与 未 飞行 的 性 成 熟 处 
女 蜂 王 mRNA 表达 差异 。 从 测序 获取 的 Tags 分 析 
结果 以 及 饱和 度 中 可 以 看 出 , 样品 的 测序 结果 较 
好 。 并 且 也 发 现 飞行 蜂王 与 未 飞行 蜂王 在 Clean 标 
签 数 、Clean 标签 种 类 以 及 相关 基因 具有 明显 的 
区 别 。 

核糖 体 和 蛋白 除 参与 核糖 体 构成 和 和 蛋 日 质 合 成 之 
外 , 还 在 复制 、 转录、RNA 加 工 等 过 程 中 发 挥 重 要 
作用 , 并 在 细胞 的 增殖 、 凋 亡 以 及 发 育 等 多 种 调控 
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方面 起 作用 ( 侯 怡 铃 等 , 2010)。 我 们 发 现 有 6 个 核 
糖 体 蛋白 基因 在 飞行 的 蜂王 中 上 调 表达 。 由 于 蜂王 
在 婚 飞 交尾 后 卵 介 被 激活 ， 卵巢 迅速 脱 大 ,此 时 蜂 
王 需 要 大 量 合成 集 晶 储存 到 卵 里 , 所 以 推测 这 些 核 
糖 体 蛋白 基因 的 上 调 表 达 , 可 能 与 此 时 蜂王 需要 大 
量 合成 蛋白 有 关 。 

细胞 色素 P450( 简称 CYP450) 是 一 类 亚 铁 血红 
素 一 一 硫 醇 盐 和 绰 日 的 超 家 族 , 它 参 与 内 源 性 物质 和 
包括 药物 、 环境 化 合 物 在 内 的 外 源 性 物质 的 代谢 。 
作为 一 种 末端 加 氧 酶 , 参与 了 生物 体内 的 省 醇 类 激 
素 合成 等 过 程 ( 李 分 和 刘 小 末 ,，2012; 首 调 义 等 ， 
2012)。 本 人 研究 发 现 3 个 CYP450 基因 在 飞行 的 蜂 
王 中 上 调 表 达 , 这 可 能 是 由 于 蜂王 为 了 吸引 雄 蜂 与 
之 交配 , 在 飞行 过 程 中 释放 更 多 的 信息 素 或 者 释放 
某 一 种 特殊 的 化 学 信息 素 , 需要 CYP450 基因 参与 
合成 。 具 体 有 待 于 进一步 研究 与 论证 。 

羽 外 , 在 这 些 差异 表达 基因 中 我 们 发 现 有 一 个 
气味 结合 生日 基因 Opp9 (CenBank 登录 号 : NM_ 
001040225. 1 ) 在 飞行 的 蜂王 中 上 调 表 达 。 气 味 结合 
蛋白 是 位 于 嗅 党 感 器 淋巴 液 内 里 的 一 种 水 溶性 的 小 
分 子 蛋 日 , 是 嗅 党 信号 传导 途径 起 始 端的 重要 功能 
组 分 。 它 们 的 功能 是 将 气味 分 子 传输 到 分 布 于 咒 觉 
神经 树 突 状 腊 上 的 感受 项 上 。 而 Opp9 基因 仅 在 较 
罕 的 发 育 历 期 中 表达 , 一 般 仅 在 蜂王 卵巢 以 及 胚胎 
早期 ( 玲 红 宽 ,2007) ,本 人 研究 发 现 0bp9 基因 在 飞 
行 蜂王 中 的 表达 量 高 于 未 飞行 的 蜂王 ,可 能 参与 蜂 
王 婚 飞 的 过 程 中 感受 周围 环境 中 气味 的 的 变化 。 

丝氨酸 和 蛋白酶 是 一 类 以 丝氨酸 为 活性 中 心 的 重 
要 蛋白 水 解 酶 , 在 生物 有 机 体 中 起 着 重要 而 广泛 的 
生理 作用 。 在 真 核 生 物 中 , 丝氨酸 重 日 酶 生理 功能 
是 间接 完成 的 , 是 通过 一 系列 的 酶 或 重 日 质 激 活 来 
完成 , 在 这 个 激活 过 程 中 , 有 些 丝 氟 酸 人 熏 日 酶 成 为 
上 游 激活 因 了 于 ， 而 另 一 些 则 成 为 下 诉 执 行 因 了 于 ( 汪 
世 华 等 , 2007) 。 本 试验 研究 发 现 , 在 性 成 熟 处 女 
蜂王 飞行 过 程 中 , 有 3 种 丝 气 酸 重 日 酶 基因 表达 上 
调 , 为 一 种 丝氨酸 入 日 酶 基因 表达 下 调 。 这 也 说 明 
这 4 种 丝氨酸 和 集 日 酶 在 蜂王 飞行 过 程 中 激活 或 抑制 
某 些 重 白 酶 原 , 具体 有 符 于 进一步 挖 气 与 验证 。 
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